Sliding mode control design for oblique wing aircraft in wing skewing process

When the wing of Oblique Wing Aircraft (OWA) is skewed, the center of gravity, inertia and aerodynamic characteristics of the aircraft all significantly change, causing an undesirable flight dynamic response, affecting the flying qualities, and even endangering the flight safety. In this study, the dynamic response of an OWA in the wing skewing process is simulated, showing that the three-axis movements of the OWA are highly coupled and present nonlinear characteristics during the wing skewing. As the roll control efficiency of the aileron decreases due to the shortened control arm in an oblique configuration, the all-moving horizontal tail is used for additional roll and the control allocation is performed based on minimum control energy. Given the properties of pitch-roll-yaw coupling and control input and state coupling, and the difficulty of establishing an accurate aerodynamic model in the wing skewing process due to unsteady aerodynamic force, a multi-loop sliding mode controller is formulated by the time-scale separation method. The closed-loop simulation results show that the asymmetric aerodynamics can be balanced and that the velocity and altitude of the aircraft maintain stable, which means that a smooth transition is obtained during the OWA’s wing skewing.     

Six sigma robust design optimization for thermal protection system of hypersonic vehicles based on successive response surface method

Lightweight design is important for the Thermal Protection System(TPS) of hypersonic vehicles in that it protects the inner structure from severe heating environment. However, due to the existence of uncertainties in material properties and geometry, it is imperative to incorporate uncertainty analysis into the design optimization to obtain reliable results. In this paper, a six sigma robust design optimization based on Successive Response Surface Method(SRSM) is established for the TPS to improve the reliability and robustness with considering the uncertainties. The uncertain parameters related to material properties and thicknesses of insulation layers are considered and characterized by random variables following normal distributions. By employing SRSM, the values of objective function and constraints are approximated by the response surfaces to reduce computational cost. The optimization is an iterative process with response surfaces updating to find the true optimal solution. The optimization of the nose cone of hypersonic vehicle cabin is provided as an example to illustrate the feasibility and effectiveness of the proposed method.     

超临界环境下煤油替代物液滴燃烧特性数值研究

为研究液氧/煤油火箭发动机燃烧室内经喷注形成的煤油液滴的燃烧过程，基于实际气体状态方程、高压热物性修正、高压气液平衡和详细化学反应动力学，建立一维的全瞬态液滴燃烧模型，对超临界环境下两组分煤油替代物液滴的燃烧特性及液滴初始直径的影响进行仿真研究。结果表明，在超临界环境下，相比于煤油液滴纯蒸发过程，煤油液滴燃烧过程的迁移时刻大大提前；煤油液滴着火之后很快进入超临界燃烧阶段，此时液滴燃烧过程可以看成中心附近的燃料高浓度区与外侧氧气高浓度区之间的扩散燃烧过程；煤油液滴的火焰半径先增大，达到最大值之后开始减小，并减小为零，火焰温度在着火之后快速上升至最大值，并基本保持不变，在火焰半径减小为零之后开始降低；随着液滴初始直径的增大，火焰特性以及液滴中心参数变化曲线趋势不变、整体延迟，着火时间、迁移时间和液滴寿命增大。     

类X-37B升力体再入飞行通信特性中断研究

针对再入非平衡等离子流场的数值仿真方法及等离子体鞘套对不同波段电磁波传输特性的预测开展研究。首先，利用某工程钝头体飞行器再入实测通信飞行测量数据，初步校验数值仿真等离子体特性和通信特性方法的正确性；其次，在典型弹道下针对升力体再入飞行器开展等离子体鞘套数值仿真，分别给出S波段、Ka波段与K波段电磁波衰减特性，研究指出结合当前国内中继星体制情况采用高频段电磁波Ka波段或K波段中继通信方式是解决通信黑障的技术途径之一；最后，研究了电磁波数衰减受传播方向角和攻角的影响规律。研究结果可为发展解决我国再入飞行器通信黑障的设计方案提供参考。     

非预混条件下的旋转爆轰燃烧室双波头演化过程数值模拟

针对旋转爆轰燃烧室双波头演化过程中流场结构变化的问题,对非预混条件下的旋转爆轰燃烧室从起爆到形成稳定的双波头过程进行了数值模拟研究。研究结果表明,从起爆到形成稳定爆轰过程,燃烧室主要经历了起爆、爆轰波对撞和稳定爆轰三个阶段；在爆轰波对撞阶段,首次对撞是两个爆轰波间的对撞,由于对撞点处缺少新鲜混合气,从而在对撞结束后衰减为两个压力波。第二次对撞是两个压力波间的对撞,因为在第二次对撞点附近存在新鲜混合气来支撑爆轰波的持续传播,故对撞结束后产生了一个爆轰波和一个较弱的压力波；第二次对撞发生后,燃烧室内的压力波反射叠加并形成局部高压区,此高压区压缩气体使气体温度升高,高温气体引燃混合气后,最终发展成为第二个爆轰波；稳定阶段,两个爆轰波均能稳定自持传播,爆轰波峰面压力可达1.45MPa,波后温度为2500K,爆轰波速度稳定在1738m/s,产生的推力与比冲分别为79.76N和2312.15s；斜激波的存在使燃烧室出口平面流场产生了较大波动。      

A comprehensive study of the 2016 Mw 6.0 Italy earthquake based on high-rate (10 Hz) GPS data

The 2016 M_w 6.0 Italy earthquake is successfully recorded by the near-field 10?Hz GPS and 200?Hz Strong Motion (SM) stations, providing valuable data for this study. A comprehensive study of this earthquake is carried out based on GPS data, which contains coseismic deformations analysis, noise analysis, seismic wave picking, and magnitude determination. The noise of most GPS-derived displacement waveforms can be described as a combination of white noise, flicker noise, and random walk noise after the earthquake occurrence, and the spectral indices vary significantly for most stations, implying that the seismic signals have affected the noise characteristic of GPS-derived displacement waveforms. S-transform is employed to assess the GPS capability to detect the seismic arrival time. The SM station AMT and the GPS station AMAT are in good agreement in seismic wave picking, and the difference is only 1.2?s in the north component, suggesting that the outcome of seismic wave picking using GPS data is reliable. Then, a classic empirical formula is employed to determine the moment magnitude. A robust moment magnitude (Mw 5.90) can be estimated by the nine GPS stations with about 23.9?s. If four GPS stations near the epicenter is chosen to determine the magnitude, it only take 13.0?s to retrieve a reliable preliminary (Mw 5.82) magnitude, which is 5.4?s ahead of nine stations. In addition, Cross Wavelet Transform (XWT) is adopted to measuring the correlation and phase relationship between GPS and SM records. The result of XWT analysis indicates 10?Hz GPS is capable of capturing reliable and accurate coseismic dynamic deformations, as evidenced by the XWT-based semblance being close to 1 between GPS and SM records. The above results confirm the capability of 10?Hz GPS to capture coseismic dynamic deformations, detect seismic arrival time, and determine earthquake magnitude. Moreover, rapid magnitude determination based on 10?Hz GPS data can be regarded as an important supplement to Earthquake Early Warning (EEW).     

预旋喷嘴径向角度对预旋特性影响的数值研究

为了降低低位进气预旋流路的气动损失，针对带有不同径向角度（0°~30°）预旋喷嘴的预旋系统进行了数值仿真，并对流动特性、温降特性和比熵增特性进行了分析。结果表明:随着预旋喷嘴径向角度的增大，预旋系统无量纲温降先增大后减小，流动阻力减小，预旋系统的流量随之增大。旋转雷诺数为2.3×107时，预旋喷嘴带径向角度的预旋系统无量纲温降比传统喷嘴最大可提高18.3%，存在某一角度使预旋温降特性达到最好。预旋系统内的耗散主要发生在预旋腔和共转盘腔内，径向角度为10°时其比熵增变化量分别占整个预旋系统总体比熵增的42.4%和30.2%；合理设计预旋喷嘴的径向角度，能改善预旋腔内气流的流动效果，并且可以减少整个预旋系统的不可逆损失。      

一种针对局部结构的有限元模型修正方法

针对实际结构有限元模型（FEM）的建模误差通常仅存在于局部区域，提出了一种对局部结构单独进行模型修正的方法。首先，根据频响函数（FRF）解耦理论得到由残余结构频响函数与包含待修正参数的局部结构动刚度所重构的整体结构频响函数的拟合值，然后通过迭代优化使其与测量值的残差最小化，从而得到参数的极大似然估计。在此基础上，将残差关于参数的灵敏度以局部结构动力学矩阵表示，建立了模型修正的基本方程，利用整体结构的测试数据即可直接对分离出来的局部结构进行模型修正。最后，对喷气式飞机和三角机翼飞机分别进行了数值模拟和实验研究，验证了所提方法的可行性和有效性。结果表明，所提方法可以成功地用于仅局部区域含有建模误差的实际结构有限元模型的修正，修正后的有限元模型的动态特性与实际结构有较好的一致性。     

面向快速响应产品的星上目标定位及识别方法

针对遥感卫星星上目标的定位及识别，本文面向快速响应应用需求，首先对星上目标定位及识别工作模式进行研究，总结出狭义和广义两种模式，其次重点对像方-物方坐标双向定位算法进行梳理，并针对典型地物目标星上自动识别算法提出研究思路，最后给出算法的应用方法流程。本文研究为快速响应产品星上智能处理提供方法参考。     

一种冲击响应谱试验设备校准装置的研制

介绍了冲击响应谱基本概念及其试验设备的分类与工作原理，参考现行冲击试验台检定规程提出了针对该设备的校准方法，设计了一种基于虚拟仪器技术的新型冲击响应谱试验设备校准装置，给出了典型冲击波形（半正弦、矩形、后峰锯齿）的冲击响应谱仿真计算结果并进行了验证，计算结果与理论结果完全一致，证明算法准确可靠，能够实现对冲击响应谱试验设备的校准。